Stability of the compressible quantum Hall state around the half-filled Landau level 
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We study the compressible states in the quantum Hall system using a mean field theory on the 
von Neumann lattice. In the lowest Landau level, a kinetic energy is generated dynamically from 
Coulomb interaction. The compressibility of the state is calculated as a function of the filling factor 
v and the width d of the spacer between the charge carrier layer and dopants. The compressibility 
becomes negative below a critical value of d and the state becomes unstable at v — 1/2. Within a 
finite range around v — 1/2, the stable compressible state exists above the critical value of d. 

PACS numbers: 73.40.Hm, 73.20.Dx 



The fractional quantum Hall effect (FQHE) jlj] is ob- 
served in a two-dimensional electron system in a strong 
magnetic field around the rational filling factor v with an 
odd denominator. The formation of the incompressible 
liquid state Q successfully explains the FQHE. Recent 
experiments of the quantum Hall system around the ra- 
tional filling factor v with an even denominator (ex. 1/2) 
indicated the formation of the Fermi surface [6™. The 
composite fermion theory |t],^| can reasonably explain 
these experiments. However, the existence of the stable 
compressible Fermi liquid at v = 1/2 is assumed a priori 
in this theory. It is more desirable to start from a micro- 
scopic theory and study the stability of the Fermi liquid 
state in the lowest Landau level space. Along this line 
a microscopic theory of the composite fermion has de- 
veloped HJTtJ an d an effort is made to obtain the finite 
compressibility at v = 1/2 p]]. Until now, understand- 
ing of the Fermi liquid property of the quantum Hall 
system in the microscopic level is unsatisfactory. The 
purpose of the present work is to study the stability of 
the compressible liquid state around v — 1/2 using a 
mean field theory in the von Neumann lattice formalism 
p2| . We show that the stable Fermi-liquid- like state is 
formed if a distance between the 2d electron layer and 
dopants is larger than a critical value. We call this state 
the quantum Hall gas (QHG) state. 

In a strong magnetic field the free energy is quenched 
and kinetic energy, pressure, and compressibility are gen- 
erated from the Coulomb interaction. Hence the pres- 
sure and compressibility of the QHG become negative in 
the jelium model where the uniform positive charge dis- 
tribution cancels the Coulomb energy between electrons 
p3[ . In realistic high-mobility samples, however, the pos- 
itive charges (dopants) are separated from the negative 
charges (carriers) by the spacer in order to obtain high- 
mobility and this separation d gives a positive contribu- 
tion to the energy of the system. Owing to this con- 
tribution, the pressure and compressibility can become 
positive for a large d and the QHG state becomes sta- 
ble. The critical values d c \ and d c i for the pressure and 
compressibility are obtained in the present paper. 



In the lowest Landau level space, the total Hamiltonian 
for a two-dimensional system of this situation is given by, 

H = ~ J Mr' : (p(r) - Po )V(r - r')(p(r') - p ) : 
+ J d 2 rd 2 r'p(r)(V(r - r') - V d (r - r'))p . (1) 

Here p(r) = ip'(r)i/j(r) and the electron field operator 
ip(r) is expanded by the Wannier basis on the von Neu- 
mann lattice HU as 

^(r) = ^6(X)(r|/ ®/3 x ), (2) 
x 

where b is the anti-commuting annihilation operator and 
X is an integer-valued two-dimensional coordinate. The 
base functions |/o <$ /?x) are orthonormal complete set 
in the lowest Landau level space. Expectation values 
of the position of the base functions are located at two- 
dimensional lattice sites a(rm, n/r) for X = (m, n) where 
a = ^j2ith/eB. The parameter r is fixed later by the 
minimum energy condition at v = 1/2. pq is the posi- 
tive background charge uniformly distributed in the sheet 
separated from the electron gas by a distance d. The 
colons in Eq. (|l]) mean that the normal ordering should 
be taken. We treat the electron as a spinless fermion. 
The potential V and Vd are given by 

V(r) = q 2 /r, V d (r) = ^/vWd 2 . (3) 

The second term of the right-hand side in Eq. (Q) gives 
2irq 2 i/ 2 d/ 'a 2 per area under the charge neutrality condi- 
tion. For simplicity we set a = 1 and h = c = 1. 

By carrying out the integration in Eq. (Q) , the Hamil- 
tonian becomes 

H = H + 2nq 2 v 2 L 2 d, (4) 

H ° = \T, : & t (X 1 )fe(Xi)nX,Y,Z)6t(X 2 )fe(X' 2 ) :, (5) 

x 

where X = X x - Xi, Y = X 2 - X 2 , Z = X x - XJ, and 
L 2 is an area of the system. The system is, then, trans- 
lationally invariant on the lattice. Therefore, if this in- 
variance is unbroken in the ground state, the formation 
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of the Fermi surface is expected. Interaction potential 
V (X, Y, Z) is given by 



d 2 k 



drpi d 2 p 2 - 



(k) exp i{pi • X 



(2n) 2 J_ w (2tt) 2 (2tt) 2 

+p 2 Y + k Z + /( Pl , k) - /(p 2 , k)}, (6) 



where u(k) = g 2 exp(-k 2 /47r)V r (k) and F(k) = 27r/|k| 
for k ^ and 1^(0) = 0, respectively. Here, k = 
(k x /r,rk y ) and Pi is the lattice momentum in the Bril- 
louin zone \pi\ < it. The phase function /(p, k) is defined 
by 

f(p, k) = -fc x (2p y + k y )/4ir - A(p + k) + A(p), (7) 

where A(p) is any real function satisfying the boundary 
condition e iA( P +27rN)-iA( P ) = (_ 1 )W,+w ! , e - l w H p,_ Her6i 

N x and N y are integers. Note that /(p, k) is writ- 
ten by using a line integral along a straight line as 
/; +fe (A(p) - 0A(p)) • dp with A(p) = ( Py /2n,0). A(p) 
represents a magnetic field in the momentum space and 
A(p) corresponds to the gauge degree of freedom. The 
Hamiltonian H is invariant under the transformation 



6(X) 



JK-Xi 



(8) 



up to the gauge transformation jL4| . Eq. (g) corresponds 
to the translation p — > p + K. This invariance is similar 
to the K- invariance in the composite fermion model [ fll| . 
In our case, on the other hand, Eq. (||) is equivalent to 
the magnetic translation in the momentum space. 

The interaction potential satisfies the following sum 
rules: 



E z nX,Y,Z) = fi(27rX)^, 

E x nx,Y-x ) z) = u (z)^ ) 



(9) 



where v(X) = jr^-t 1 / 2 '^ / ((tt/2)X 2 ), which is the 
Fourier transform of v, where Iq is the modified Bessel 
function of the first kind and X = (rm,n/r) for X = 
(m,n). X is a position of the Wannier basis in the real 
space. 

We introduce a mean field for the QHG state which has 
the translational invariance on the von Neumann lattice. 
Namely, 



[/ (X-X) = (fe t (X)6(X)), 



(10) 



and U (0) = v. If U (X) ^ for X ^ 0, the symmetry 
under Eq. (||) is broken. Using the sum rules of Eq. (||), 
the mean field Hamiltonian for Ho becomes 

H mf = J2 {-^o(X - X>(X' - X)fet( X )6(X) (11) 

X,X' 

+^U (X - X')w(X' - X)(7 (X' - X)}. 



The first term of the right-hand side is a kinetic energy 
dynamically induced from the Coulomb interaction. It 
is remarkable that the continuum system with a mag- 
netic field is transformed to the lattice system without 
a magnetic field. In the composite fermion model, the 
external magnetic field is canceled by the Chern-Simons 
gauge field at v — 1/2. In our model, on the other hand, 
the magnetic field in the momentum space disappears for 
a mean field Uq(X. — X') at arbitrary filling factor. 

By Fourier transforming Eq. (pl]), and using Eq. (|lC|), 
we obtain self-consistency equations for the one-particle 
kinetic energy e(p) in the Hartree-Fock approximation as 



( \ f d2 Pi - 



vt (Pi - p)0{na{v) ~ e{px,v)), 



dp 



e(no(v) - e(p,i>)), 



(12) 



where /io is the chemical potential for the mean field 
and V{ is defined by "Df(p) = X)n^(p + 27rN). Equa- 
tion (|l^) determines a self-consistent Fermi surface. Ex- 
istence of the Fermi surface breaks K-invariance. We ob- 
tain the solution that breaks if^-invariance but preserves 
if^-invariance. We use the lines \p y \ — nv as the self- 
consistent Fermi surface which is a solution of Eq. (jT|) 
for an arbitrary v | fl5| . The one-particle states are filled 
in \p y \ < 7tza Then the kinetic energy e(p, v) becomes 
constant in the p^-direction and is given by 

e(P,^) =s(p y ,v) =-J -^vxipy-p), (13) 
Viip) = q 2 r e- r2{p+2m)2/ ^K {r 2 {p + 27m) 2 /8^), 

n 

where Kq is the modified Bessel function of the second 
kind and the chemical potential /io is written as {y) = 
£(kv,v). The total chemical potential fi is calculated by 
the total energy per particle E(y) as /i = E(v) + vE'iv). 
The kinetic energy e{p yi v) is obtained numerically and 
is shown in Fig. 1 for v = 1/4, 1/2, and 3/4. 

Owing to the Coulomb interaction, the Fermi velocity 
de/dpy diverges logarithmically at the Fermi momentum 
p y = ±7riA Hence, there exists the Coulomb gap in the 
one-particle spectrum. This is consistent with the exper- 
iment |l(| which suggests that there is no true metallic 
phase around v = 1/2. Naively it seems that the diver- 
gence of the Fermi velocity makes the compressibility n 
vanishes because k is proportional to the inverse of the 
Fermi velocity in the ordinary electron gas. As is seen 
later, however, the v dependence of e{p yi v) contributes 
to k negatively. Furthermore, the distant background 
contribution to k is positive. Then these two contribu- 
tions compete with each other and cause a transition. 

First, we calculate the total energy of the QHG state. 
The energy per particle due to mean field Hamiltonian 
H mf is given by 



2 



£mfH = ^ ^e{p y ,u). (14) 

We calculate the fluctuation around the mean field and 
find the Hartree term contributes to the energy per par- 
ticle in the lowest order as 

E h (v) = ^Y. EMXM2vrX)[/o(-X) (15) 
x 

= rL 2^F^ — ■ 

n=l 

We approximate the total energy per particle E(y) as 

E{v) = E (v) + 2irq 2 vd, (16) 
E Q {v)=E ai {y)+E b .{v). 

Eq(1/2) becomes minimum at r = 1.64, which is used in 
this letter. The minimum value of Eq{1/2) is — 1.086g 2 
(= -0.433<2 2 /'b, l B = y/h/cB). This value is slightly 
higher than values of the charge density wave (CDW) 
calculation |l7]] which is shown in Fig. 2 by a dot. In 
the CDW, the higher order correction is small because of 
the energy gap. In the QHG, however, the higher order 
correction might be large compared with the CDW. The 
value of Eh is positive and is smaller than one tenth of 
the magnitude of E m {. Although E^ has a small contri- 
bution to the energy, this term is necessary for £^(1/2) 
to have the minimum as a function of r because E m { and 
£"h are monotonically decreasing and increasing function 
of r, respectively, at fixed filling factor. In Fig. 2, E{v) 
is plotted for d — 0, 0.1, and 0.2. At v = 1/3, the value 
of Laughlin's incompressible state is shown in Fig. 2 by 
a cross. The energy of the QHG for d = is substan- 
tially higher than the Laughlin's value. Therefore the 
incompressible state is preferable to the QHG state. 

Second, we calculate the pressure of the QHG state. 
The pressure P is defined by the derivative of the to- 
tal energy per particle as P = v 2 E'(v). In Fig. 3 the 
pressure is plotted for d = 0, 0.1, and 0.2. Using the 
particle-hole transformation b, it can be shown that 
£'(1/2) obeys the relation £'(1/2) = -2(£(l/2)-£(l)), 
where E(l) = —nq 2 /2 + 2nq 2 d. Therefore the pres- 
sure at the half-filling becomes negative below a critical 
value d = d c i, which is determined by £(1/2) = E(l) as 
del = (£b(l/2) + irq 2 /2)/nq 2 = 0.154 in the unit of a. 
For v = 1/4 and 3/4, the critical values are numerically 
obtained as 0.253 and 0.138 respectively. As is seen in 
Fig. 2, E(y) is a decreasing function for d = 0. This 
behavior was verified by the exact diagonalization of the 
small system [|l8|. Thus, if E{v) is a smooth function, 
the pressure is negative for d = 0. The negative pressure 
does not always mean the instability of the gas state. For 
the stability, the positive compressibility is necessary. 

Finally, let us calculate the compressibility of the QHG 
state. The compressibility k is defined by the derivative 



of the pressure P{v) as 1/k = vP'{v). In calculating k 
it is rather convenient to use the relation 1/k = v 2 ijl{y). 
Using Eq. @ we obtain 

00 

H'(v) = -Vi{2ttv) + 2q 2 re-^ n/r)2 cos{2-KVn)/n 

n=l 

+Anq 2 d. (17) 

The first term of the r.h.s. comes from /j,q and the sec- 
ond term comes from E^. In Fig. 4, l//j/(y), which is 
proportional to k, is plotted for d = 0, 0.1, and 0.2. The 
positive k region appears symmetrically with respect to 
v = 1/2 for d > d c2 . The critical value d c2 at v = 1/2 
is determined by //(1/2) = as d c2 = 0.115 in the unit 
of a. For v = 1/4 and 3/4, the critical values are also 
obtained as 0.151. 

The negative compressibility has been observed by us- 
ing the electric field penetrating technique |fl9|| . This 
method makes it possible to measure the compressibility 
of the two-dimensional electron gas discriminating the 
background charge contribution. The observed negative 
compressibility may correspond to the Eq. ( |l7| ) for d = 0. 
Near v = and 1 the compressibility is always nega- 
tive in Fig. 4. A transition between the QHG and the 
Wigner crystal is expected when the k diverges. Note 
that d c i > d C 2 at v — 1/4 and 1/2, and d c \ < d c2 at 
v = 3/4. At B ps 10 T, d c i and d c2 are on the order of 
1 nm. By varying B, d c \ and d c2 vary as oc B~ x / 2 . In 
the high- mobility GaAs-AlGaAs heterostructure, d >10 
nm, the QHG state becomes stable. Changing a mag- 
netic field or a width of the spacer, the transition could 
be observed in experiments. The QHG state might be 
changed into the incompressible liquid state or the insu- 
lator phase. 

In the composite fermion model, the Fermi wave num- 
ber k,F at v = 1/2 is given by kp = y/<iirn e , n e is the 
density of the electrons, which was verified by experi- 
ments [|]-[|. In our model, the Fermi wave number pp 
on the von Neumann lattice is given by pp — nv, which 
corresponds to the the Fermi wave number kp on the real 
space given by kp = m/r/a — my/Un^. For r = 1.64 and 
v = 1/2, rix^/v — 3.64, which is very close to the value 
V4^ = 3.54. 

In summary, we study the QHG state at arbitrary v us- 
ing the von Neumann lattice. The QHG state has a Fermi 
surface and is stabilized by the spacer around v = 1/2 
for d > d c2 and a transition occurs for d < d c2 in the 
quantum Hall system without disorder. 
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Fig 1. Kinetic energy e{p y ,v) for v — 1/4, 1/2, and 
3/4. The unit of the energy is q 2 /a. The gradient of e 
diverges at p y = ±ttza 
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Fig. 4. Compressibility k{v) times v 2 for d = (dashed 
line), 0.1 (dash-dotted line), and 0.2 (solid line). The unit 
of the k is a 3 /q 2 . 
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Fig. 2. Energy per particle E(v) for d = (dashed 
line), 0.1 ( dash-dotted line), and 0.2 (solid line). The 
unit of the energy is q 2 /a. The cross (x) shows Laugh- 
lin's value at v = 1/3 and dot (•) shows CDW value at 
^ = 1/2. 
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Fig 3. Pressure P{v) for d — (dashed line), 0.1 (dash- 
dotted line), and 0.2 (solid line). The unit of the pressure 
is q 2 /a 3 . 
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